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(N ■ ABSTRACT 

o : 
o . 

I Narrow-band searches for Lyman-a emission are an efficient way of identifying star- 
forming galaxies at high redshifts. We present Keck telescope spectra confirming red- 

■ shifts z ~ 5.7 for three objects discovered in the Large Area Lyman Alpha (LALA) 
If^ • survey at Kitt Peak National Observatory. All three spectra show strong, narrow emis- 

. sion lines with the asymmetric profile that is characteristically produced in high redshift 

y—i . Lyman-a emitters by preferential HI absorption in the blue wing of the line. These ob- 

^ i jects are undetected in deep B]y, V, R, and A ~ 66OOA narrow-band images from the 

Ti^lj- I NOAO Deep Wide-Field Survey and from LALA, as expected from Lyman break and 

I Lyman-a forest absorption at redshift z ~ 5.7. All three objects show large equiv- 

' alent widths (> 150A in the rest-frame), suggesting at least one of the following: a 

■ top-heavy initial mass function, very low stellar metallicity, or the presence of an active 

■ nucleus. We consider the case for an active nucleus to be weak in all three objects due 
Q^i to the limited width of the Lyman-a emission line (< 500km s^^) and the absence of 
Q ! any other indicator of quasar activity. The three confirmed high redshift objects were 
|H I among four spectroscopically observed targets drawn from the sample of 18 candidates 
^ I presented by Rhoads & Malhotra (2001). Thus, these spectra support the Lyman-a 

' emitter population statistics from our earlier photometric study, which imply little evo- 
lution in number density from z = 5.7 to z = 4.5 and provide strong evidence that the 
reionization redshift z,. > 5.7. 
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1. Introduction 

Narrow-band searches were originally proposed as a method of finding young, high redshift 
galaxies in their first throes of star-formation (Partridge &; Peebles 1967). This method has now 
proven its worth in discovering star-forming galaxies at ever increasing redshifts. A large fraction 
(9 out of 11) of known z > 5 galaxies have strong Lyman-a emission (Dey et al. 1998, Dawson 
et al. 2002, Weyman et al. 1998, Hu ct al. 1999, Hu ct al. 2002, this paper), the only exception 
being a pair of Lyman break galaxies at z=5.34 (Spinrad et al 1998). This is partly due to the 
ease of confirming galaxies with strong emission lines and partly because one can avoid night sky 
lines by using narrow-band selection methods. Ultimately, the Lyman-a emitters at high redshifts 
will provide useful constraints on models of galaxy formation and evolution. A more immediate 
application is to constrain the epoch of reionization. In a neutral universe, Lyman-a photons from 
galaxies are resonantly scattered into large, diffuse halos, whose extremely low surface brightness 
will make their detection very challenging. Thus, the transition to a neutral universe will be marked 
by a sharp drop in Lyman-a galaxy counts. 

In the Large Area Lyman Alpha (LALA) survey (Rhoads et al. 2000, Paper I) wc use narrow- 
band imaging with the wide-field Mosaic cameras at Kitt Peak National Observatory (KPNO) and 
Cerro Tololo Inter-American Observatory (CTIO) 4-meter telescopes to image a relatively large 
area (36' x 36') of the sky at once. This allows us to obtain good statistics on the candidates 
selected by imaging, and use them to study galaxy formation and evolution at high redshifts. 
Confirming candidates from narrow-band imaging by spectroscopy requires large telescopes and is 
a slower process. But as we show in this and companion papers, the success rate for spectroscopic 
confirmations is high, so conclusions drawn from mere imaging candidates are probably accurate 
(Rhoads & Malhotra 2001, hereafter RMOl; Malhotra & Rhoads 2002; Ouchi et al. 2002). 

2. Observations 

2.1. Narrow- band and Broad-band Imaging 

The LALA survey concentrates on two primary fields in constellations Bootes (J142557-I-3532) 
and Cetus (J020520-0455). Each field is 36 x 36 square-arcminutes in size, corresponding to a 
single pointing of the 8192x8192 pixel NOAO Mosaic CCD cameras that we have used for all of 
our imaging. The LALA fields are located within the NOAO Deep Wide-Field Survey (Jannuzi 
Sz Dey 1999), which provides deep broad-band images in Bw, R, and I filters that support our 
narrow-band data. 

The z 5.7 search has been conducted in the Bootes field using two narrow-band filters with 

central wavelengths Ac of 815 and 823 nm, full widths at half maximum (FWHM) transmission 
7.5nm, and peak throughputs ~ 90%. The wavelengths were chosen to fall in a gap in the night sky 
emission line spectrum, thus minimizing sky noise (Figure 1). The filter bandpasses deteriorated 
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below these specifications outside a central circle of ~ 30' diameter, resulting in an effective search 
area of ~ 710 square arcminutes. Within this area, we identified a sample of 18 candidate Lyman-a 
emitters based on strong (> 5a) detections in a narrow-band filter that substantially exceed the 
continuum flux level (observer- frame equivalent width > 80A and a 4c7 excess over / band), and 
non-detections (< 2a) in NDWFS B^r and LALA V filter images. This sample is presented in 
RMOl, to which we refer the reader for more details on the photometry and selection criteria. 

2.2. Spectroscopic Observations 

Spectroscopic followup of four LALA z = 5.7 candidates was conducted at the Keck I 10m 
telescope on UT 2002 March 12, 13, and 14. The observations were made using the Low Resolution 
Imaging Spectrometer (LRIS; Oke et al. 1995) equipped with the 400 lines/mm grating (A5 = 
8500 A) and multislit masks with 1'.'4 wide slitlets. The 4 LALA z=5.7 candidates were targeted 
using 2 multislit masks; LALA5 1-05 and 1-06 were observed on March 12 for ~ 110 minutes, 
and LALA5 1-01 and 1-03 were observed on March 13 and 14 for a total of 175 minutes. All 
observations were made with the slitlets oriented at roughly 320° (east of north). Seeing on the 
nights varied between 0^7 to I'.'O. Because the slits were wider than the seeing FWHM, spectral 
resolution becomes a function of source size, from ^ 6.5A for unresolved sources in the best seeing 
data to I2A for sources that fill the slit. The LALA targets on these masks arc all < 1" in size 
(based on the narrow-band detection images) implying an effective resolution in the range 6 to 9A. 

These data were reduced using a combination of standard IRAF ONEDSPEC and TWOD- 
SPEC tasks together with homegrown slitmask reduction tasks developed by Stern, Bunker, & 
Stanford (private communication) for reducing LRIS slit mask data.^ 

3. Results 

Three of the four z ^ 5.7 candidates were confirmed as Lyman-a emitters with the expected 
redshifts and line fluxes. We designate these as LALA J142546.76-F352036.3, LALA J142647.16-F353612.6, 

and LALA J142630.34-F354022.5. These correspond respectively to candidate IDs LALA5 1-03, 
LALA5 1-06, and LALA5 1-05 from RMOl. The extracted spectra of the candidates are presented 
in Figure 1. The photometric and spectroscopic properties of these objects are summarized in 
tables 1 and 2. 

The fourth candidate on the slit masks (LALA5 1-01) had the lowest narrow-band flux among 
the four, corresponding to a line flux of (1.57 ± 0.27) x 10~^^ erg cm^^ s^^. This candidate was not 
detected in the spectra to a 5a upper limit of ~ 0.7 x 10""^^ erg cm~^ s~^. Since it was a single-filter 
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detection in the photometric survey, it was most hkely a noise spike in the narrow-band image 
and not a real object. Given the large number (~ lO^"'') of independent resolution elements in the 
images, we expect about 1 false positive at the 5cj level per LALA imaging field per narrow-band 
filter, regardless of the number of true emission line galaxies present in the survey. The number 
of false positives could be larger if the noise properties of the image are not precisely gaussian. 
A transient origin for this source is unlikely, as it is not detected in the z' filter image that was 
obtained during the same 4 night observing run as the 815 and 823 nm images. 

All of the confirmed sources show asymmetric spectral lines, with a sharp drop on the blue 
side and a more extended wing on the red side. This asymmetry is characteristic of high redshift 
Lyman-a and is physically associated with absorption by neutral hydrogen in outflows from actively 
star-forming galaxies (e.g., Stern & Spinrad 1999; Manning et al 2000; Dawson et al 2002). This 
asymmetry is good evidence against an [O II] AA3726,3729 line identification, because the close 
doublet structure of [O II] AA3726,3729 is expected to give the opposite asymmetry from Lyman-a 
under normal physical conditions for [O II] . None of the sources show any clearly detected spectral 
line besides the primary A 8200A line. (The coverage is 6850A < A < 10640A for LALA5 1-05 
and 1-06, and 5584A < A < 9400 A for LALA5 1-03.) This immediately rules out the [O III] 
AA4959,5007 doublet as a possible contaminant. 

Because the slit size exceeds the apparent angular size of the sources, the measured line widths 
contain contributions from the physical line width, the angular size of the source, and the seeing. 
An upper limit to the physical line width can be obtained by assuming the objects are unresolved 
at the 0^'7 seeing of the spectroscopic observations, providing an effective instrumental resolution 
of ^ 6A. The resulting line widths are given in Table 2. All are < 500 km s~^. 

None of the three sources shows detectable continuum on the blue side of the line in the spectra. 
There are marginal detections of the continuum emission on the red side of the line in two of the 
three sources in our spectra, with formal significance levels of 2a and 4(7. Thus, in two of the three 
sources we have evidence for a continuum break across the line, again supporting the identification 
as Lyman-a. These continuum measurements were done using a variance- weighted average of the 
flux in 1200A wide windows on either side of the lines. These quoted significance levels are based 
only on random (photon counting) error, and should be regarded as upper bounds on the true 
significance, since systematic errors will increase the uncertainty in the continuum flux. 

4. Discussion 

We now consider the physical implications of our observations, beginning with the properties 
of the individual sources and then moving on to properties of the galaxy population and of the 
universe at 2; 5.7. 

The three spectroscopically confirmed Lyman-a sources presented here all show very large 
equivalent widths, with 2a lower bounds of 150A (rest frame) and up. Such large equivalent widths 
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axe very hard to achieve with only the ionizing flux from a normal stellar population, which can 
provide an equivalent width beyond 200A for only a few xlO^ years (Chariot & Fall 1993). To 
produce equivalent widths beyond ~ 200A, the ultraviolet radiation in the galaxies must be harder 
than usual. This can be achieved in galaxies with active nuclei, with a stellar metallicity near zero 
{Z < 0.05Zq, Tumlinson k ShuU 2000; Bromm, Kudritzki & Loeb 2001), or with a stellar initial 
mass function dominated by high mass stars (Malhotra k, Rhoads 2002; Kudritzki et al. 2000). 
However, conventional (broad-lined) active galactic nuclei are excluded by the observed Lyman-a 
line widths in these sources, which have upper bounds from 220 to 460km s^^ (Table 2). Thus, if 
the sources are AGN, they must be type II AGNs, i.e., systems whose broad-line regions are hidden 
from direct view. Even type II AGN generally have larger Lyman-ct line widths than we observe. 
For radio galaxies, a composite spectrum (Stern et al 1999) shows FWHM = 1130km s~^, while 
an individual case (4c41.17; Dey et al 1997) shows a narrow component of 613 km s~^ but also a 
broad component of similar line flux and FWHM = 1373 km s~^. For X-ray selected type II AGN, 
we have one measurement at 1520 km s~^ (Stern et al 2002) and one upper limit at < 1130 km s~^ 
(Norman et al 2002). Thus, line widths seem to favor star formation activity over active nuclei as 
the driver for the Lyman-a line in these LALA sources. 

The remaining explanations for the observed equivalent widths both suggest star formation 
in primitive conditions. Low metallicity stars obviously require low metallicity gas, while top- 
heavy initial mass functions may be a natural consequence of the way low-metallicity gas {Z/Zq < 
5 X 10~^) cools and fragments (Bromm et al 2001). A final possibility is that anisotropy in the 
Lyman-a radiation biases the equivalent width measurements substantially (see Chariot and Fall 
1993; Neufeld 1991), but this explanation would require a very special source geometry. Further 
tests of the narrow-lined AGN hypothesis are possible, including searches for X-ray emission and 
for other optical emission lines characteristic of AGN (especially CIV and Hell). Such tests are 
underway for the LALA survey's larger z ^ 4.5 sample. To date we have seen no obvious evidence 
for AGN activity in any of our confirmed high-redshift LALA spectra. 

The confirmation of 3 z 5.7 Lyman-a emitters among 4 candidates observed allows us to 
refine the statistics of this population that were presented in RMOl. There we presented a sample 
of 18 candidates, based on detection significance, emission line strength, and an absence of blue 
fiux. A complication in candidate selection is that the measured equivalent widths of high redshift 
objects may be changed by intergalactic hydrogen absorption of either continuum or line flux. A 
"worst case" (and somewhat unrealistic) scenario is where the broad-band flux is substantially 
reduced while the narrow-band flux is unaffected. If we assume this case and correct the / band 
fluxes of our 18 candidates for intergalactic hydrogen absorption using a factor of 2.1 (based on the 
prescription of Madau 1995), we flnd that 13 candidates still pass the equivalent width cutoff for 
selection. The remaining 5 candidates could in principle be continuum sources at redshift just below 
5.7, where an apparent narrowband excess of a factor ~ 2 could be produced by IGM hydrogen 
absorption in the blue side of the broad band filter. We will call the 18 candidates the "full sample" 
and the 13 remaining after our worst case IGM absorption correction the "reduced sample" in the 
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foUowing discussion. All four candidates with Keck spectra belong to the reduced sample. 

Applying a 3/4 correction factor to the observed source counts for the spectroscopic success 
rate presented here, we find a number density of 540 (390) Lyman-a emitters per square degree 
per unit redshift for the full (reduced) sample. These counts are above a detection threshold of 
(1.47 ± 0.08) X 10~^^ ergcm"^ s~^, where the error bar indicates the range in detection threshold 
between the two narrow-band filters. 

The implied Lyman-a luminosity density at z ^ 5.7 then becomes 1.26 X 10^ (0.80 X 10^) 
Lo/Mpc^ for sources with Li^ya > 1-36 x IO^Lq, based on the full (reduced) sample. (Here and 
throughout this section we assume a cosmology with Hq = 70 kms~^ Mpc~^, Qm = 0.3, and 
Cl\ = 0.7.) This number is a lower bound, and the total luminosity density could be substantially 
increased if the luminosity function rises steeply below our detection threshold. 

These Lyman-a luminosities can be translated into star formation rates under an assumed 
model for the stellar initial mass function (IMF) and for the escape fraction of Lyman-a photons 
from their parent galaxies. To enable comparison with other work, we follow Hu, McMahon, & 
Cowie (1999) and apply a conversion factor of IMq/jt = lO^^ergs""^. This factor is obtained by 
taking the Lyman-a to Ha ratio for Case B recombination and no dust, together with the Kennicutt 
(1983) conversion between Ha luminosity and star formation rate. This corresponds to an IMF 
slope 2.5 (close to the Salpeter value of 2.35) over the range 1 < M/Mq < 100, continued to a 
lower mass cutoff of O.IMq with slope 1.4, with solar metallicity for all stars. We hasten to add 
that these assumptions appear inconsistent with the observed equivalent width distribution of the 
LALA Lyman-a emitters (Malhotra & Rhoads 2002), but as we cannot uniquely determine the 
correct IMF and metallicity from the available data, it is not clear how wc should replace these 
assumptions. The resulting star formation rate density aX z 5.7 is 5(3) x lO~^M0/yr/Mpc^ in 
objects with LLya > 7x 10^^ ergs~^, for the full (reduced) samples after correction for spectroscopic 
completeness. The star formation rates per object range from 7 to 14M0/yr. 

The Lyman-a luminosities allow us to constrain the metal production rate density in these 
objects much more reliably than the total star formation rate density. The ratio of ionizing photon 

production to metal production depends only weakly on the details of a stellar population, because 
the most massive stars dominate both production rates (e.g., Madau & Shull 1996). Every observed 
Lyman-a photon implies the production of > 1.5 ionizing photons (where equality is expected in the 
limit of zero dust attenuation) . The metal production rate is approximately dMz /dt = 500Li / , 
where Lj is the luminosity in ionizing photons and Mz is the mass in elements with atomic number 
Z > 6 (Madau &; Shull 1996). Thus, our Lyman-a luminosity density implies a metal production 
rate density > 8.6(5.5)M0/Myr/Mpc^ for the full (reduced) sample. Malhotra & Rhoads (2002) 
suggest that the Lyman-a emitters in the LALA sample suffer little dust attenuation, based on 
their large observed equivalent widths. Our lower bound is then likely to be near the total metal 
production in Lyman-a emitting galaxies. However, Malhotra &; Rhoads (2002) also note that 
additional objects (accounting for ^ 90% of all galaxies) could be removed from the sample by dust 
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obscuration. This would result in a total metal production rate ^ 10 times larger. For comparison, 

we estimate the metal production rate density in Lyman break galaxies at z ~ 4 to be ~ 14 times 
larger than that in our observed Lyman-a sample, based on the UV luminosity densities in Steidel 
et al. (1999). Earlier work by Madau et al (1996), adjusted to our assumed cosmology, gives 
85M0/Myr/Mpc^ at z ^ 4, for continuum selected 450 nm dropouts in the Hubble Deep Field. 
This is again a lower bound that can be increased by corrections for low luminosity sources and for 
dust absorption. 

Our spectroscopic results also reinforce the conclusion of RMOl that the space number density 
of Lyman-Q! emitters with line luminosity ivLya ^ 7 x 10^^ erg s"-*^ shows little evolution between 
2; ~ 4.5 and z ~ 5.7. In contrast, the surface density per unit redshift of Lyman-a emitters above 
observed line flux ~ 1.5 x 10^^^ ergcm"^ decreases by a factor of > 4 over the same range. This 
is in contrast to the baseline model of Haiman &: Spaans (1999), which would predict a decrease 
by at most a factor of 2 over this redshift range. 

Finally, the spectroscopic confirmation of these low- luminosity Lyman-a sources at z « 5.7 
secures the conclusion (RMOl) that the reionization redshift Zr > 5.7. The basic physics underlying 
this conclusion is simple. Lyman-a photons propagating through intergalactic gas with a substantial 
neutral fraction are subject to resonant scattering by atomic hydrogen. Thus, even flux on the red 
side of line center can be hidden by the red damping wing of a substantially neutral IGM (Miralda- 
Escude 1998; Miralda-Escudc & Recs 1998). This would obscure the Lyman-a line from view, so 
that detection of the Lyman-a line at some redshift z implies a predominantly ionized medium 
around the line source, and suggests Zr > z (Haiman &; Spaans 1999). While the resonantly 
scattered photons will eventually redshift in the Hubble flow and escape to infinity, their effective 
photosphere (or last scattering surface) is expected to subtend an angle ~ 15" (Loeb & Rybicki 
1999), reducing their surface brightness by a factor > 100 and thereby rendering them undetectable 
to surveys like LALA. 

This test was first applied by RMOl to show that Zr > 5.7, and subsequently by Hu et al. 
(2002) to argue for Zr > 6.56. To make these reionization limits firm, several nuances in the above 
physical argument deserve close consideration. In particular, these objects have to produce ionizing 
photons in order to generate the Lyman-a line by recombination. If the escape fraction /esc for 
the ionizing flux is large, the source will be surrounded by an ionized bubble in the IGM. The size 
of this bubble is limited by the (time-integrated) ionizing photon production of the source. RMOl 
show that the ionized bubble (in the dust-free case) will be large enough for Lyman-a photons 
to escape freely (i.e., 1.2Mpc) if ^43*8/080/(1 — /esc) ^ 5, where L^y^a = 10'^^L43 erg s^-*^ and the 
source is lO^ts years old, but that i43*8/esc/(l ^ /esc) ^ 2 for the observed z 5.7 LALA sources. 
The critical constraint on ig and /esc is derived from the equivalent widths of the Lyman-a sources, 
which can only achieve their large observed values in young galaxies with modest escape fractions. 
The Hu et al. (2002) source, HCM 6A, has a more modest equivalent width (25A rest frame) that 
could be produced in systems with tg ~ 1 and /esc ~ 0.75, but also shows lower line luminosity 
-L43 ~ 0.3. It therefore fulfills L^^tsfesc/ (l~/esc) ^ 8, where the last inequality comes from the 
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age of the universe at this redshift. This should be compared with a bound i^43t8/esc/(l — /esc) ^ 7 
to produce a Stromgren sphere of radius 1.2Mpc in the denser IGM at z = 6.56. Thus, it again 
appears unlikely that HCM 6A could ionize a sufficiently large bubble in a fully neutral IGM to 
allow its Lyman-a photons to escape freely to infinity. It is possible, of course, for any particular 
Lyman-a source near Zr to be more luminous than we think and have its Lyman-a flux substantially 
attenuated an effective Lyman-a optical depth of a few. Thus, HCM 6A al object may be consistent 
with Zr ^ 6 (e.g. Haiman 2002). However, the Lyman-a luminosity function is quite steep at these 
flTix levels, so the number of detectable Lyman-a sources is still expected to drop sharply at Zr. 
Haiman (2002) estimates that the suppression of observed Lyman-a flux would be a factor 10 if 
HCM 6a is in a fully neutral IGM. This would imply an intrinsic line luminosity ~ 3 x 10^"^ erg s^^, 
making HCM 6a more luminous than any object found by the LALA survey in a volume ~ 30x 
greater than surveyed by Hu et al. We therefore regard this possibility as unlikely. 

Other possible loopholes involve the detailed kinematics of the Lyman-a emitters and the 
neutral IGM, including the possible effects of gas infall, galactic peculiar motions, and galactic 
winds. Detailed consideration of these effects is beyond the scope of the current paper. However, 
we can make quick estimates of these effects. First, infall of the IGM to a newly forming galaxy 
can blueshift the gas nearest the galaxy outside the Lyman-a line core, but there must be some 
point along the line of sight where infalling gas gives way to the Hubble flow and the IGM has 
zero velocity relative to the emitted line photons. Hence, infall will not greatly reduce the resonant 
scattering opacity. Next, note that the required velocity scale to allow free escape of Lyman- 
a is 1000 km s~^. Peculiar velocities of this magnitude are restricted to the richest clusters 
of galaxies today, and should be vanishingly rare in the z ~ 6 universe. Finally, a spherically 
symmetric galactic wind with (peculiar) velocity 1000 km s~^ at the source, smoothly decreasing 
to match the Hubble flow with zero peculiar velocity at ~ 1.2Mpc, would require a total energy 
~ lQ6i-5gj.g gpj-ead over a Hubble time at 2; ~ 6, this corresponds to a mechanical luminosity of 
~ substantially exceeding the total UV/optical light output of these objects (< IO^^Lq). 

The required energy would of course be reduced in an IGM with substantial clumping or ionization. 
However, mechanical energy is typically 1% and at most ~ 10% of the starlight in a starburst 
(Leitherer & Heckman 1995), a factor of > 100 below the naive calculation above. Thus, while 
galactic winds are common features of starbursts and probably play a role in Lyman-a escape 
from local galaxies (Kunth et al 1998), they may not be able to "save" Lyman-a emission in a 
substantially neutral IGM. 

A final concern is the possibility that a low-luminosity Lyman-a source may be visible at z > Zr 
if it is located within the ionized bubble of a more luminous neighbor. In this case, the neighbor 

might be visible to the Lyman-a survey itself or to a Lyman break survey of the region (such as the 
NDWFS). The required luminosity corresponds to Iab ^ 24 — 7.5 * lg{R/2Mpc), for a Stromgren 
sphere radius R and a normal starburst ionizing spectrum. This number is only a lower bound 
because it assumes /esc = Ij which is almost certainly false for the ionizing photons. A search for 
NDWFS sources with Iab < 24 and nondetections at the 2a level in both R and Bw filters finds no 
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viable candidates within 2Mpc (projected physical separation) of any of the three spectroscopically 
confirmed sources. 

In any case, Stromgren spheres of bright neighbors may explain an isolated Lyman-a emitter 
prior to reionization (such as the Hu ct al. source at z = 6.56), but these Stromgren spheres include 
a minority of the volume before the overlap phase of reionization. Thus, a statistical sample of 
Lyman-a emitters such as the LALA survey provides should still see a drop in counts by a factor 
of > 2 at the overlap phase. The fact that we see no such drop in our z = 5.7 sample (RMOl) 
reinforces the argument that Zr > 5.7 on the line of sight to the LALA Bootes field. 

The neutral fraction required to effectively scatter most of the Lyman-a flux from low- luminosity 
Lyman-a sources can be estimated by setting r 1 at a velocity lOOAuioo km s~^ on the red 
side of line center. RMOl find that this requires a homogeneously mixed neutral fraction 
0.lAfioo[(l + z) /6.7]~^^'^ . Figure 2 summarizes our present observational understanding of reion- 
ization by combining this estimate for the neutral fraction limits placed by Lyman-a sources with 
the mass-averaged neutral fraction evolution that Fan et al (2002) infer from Gunn- Peterson trough 
searches in high redshift quasars (Becker et al 2001; Pentericci et al 2002). More detailed mod- 
eling that combines the Lyman-a luminosity function with a more realistic, inhomogeneous IGM 
is planned and will reduce the uncertainties in figure 2. Still, because the Gunn-Peterson trough 
cannot effectively probe the xhi ^ 0.01 regime, it is already clear that the Lyman-a emitting 
galaxies have a role to play in extending our knowledge of reionization to z > 6.5. 

To summarize, we have presented spectra confirming redshifts z « 5.7 for three Lyman-a 
emitting galaxies selected from the Large Area Lyman Alpha (LALA) survey. These objects were 
among four candidates observed from a total sample of 18 candidates. This implies a total sample 
of ~ 13 real z f« 5.7 galaxies, easily the largest uniformly selected galaxy sample at z > 5.5. These 
Lyman-a emitters all have narrow lines (inconsistent with broad-lined quasars) and very weak to 
absent continuum emission. They belong to the class of high equivalent width Lyman-a sources 
found in the LALA survey at z ~ 4.5. The population properties of these high equivalent width 
objects are reasonably explained by young galaxies with comparatively low metal and dust content 
(Malhotra & Rhoads 2002). The z « 5.7 population is similar to the z « 4.5 one in its number 
counts. This suggests that the sources reside in a predominantly ionized universe and hence that 
the reionization redshift Zr > 5.7. 
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Table 1. 


Photometry 






ID 


Photometric W5^^'*(2cr) 
line flux'' limit (A) 


I flux 


z' flux 

(/xJy) 


LALA J142546.76+352036.3 

LALA J142647.16+353612.6 
LALA J142630.34+354022.5 


L82 ± 0.27 
3.67 ±0.31 
2.32 ± 0.32 


> 1650 

> 195 

> 150 


-0.07 ± 0.05 
0.01 ±0.05 
-0.01 ± 0.04 


-0.10 ±0.11 
0.34 ±0.11 
0.16 ±0.11 



^Units of 10 erg cm ^. 

Note. — Photometric properties of the three confirmed z ~ 5.7 LALA sources. The 2a" 
lower bounds on equivalent width are derived from the narrow- and broad-band photom- 
etry and corrected to the spectroscopically determined redshift. They are not corrected 
for IGM absorption, which will affect both the broad and narrow-band photometry ap- 
proximately equally if the unabsorbed line profiles are centered on zero velocity. 
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Table 2. Spectroscopy 



ID 


Redshift 


Spectroscopic 
line flux^ 


FWHM 

(A) 


Av 
(km/s) 


Continuum (/iJy)*^ 
Blue side Red side 


LALA J142546.76+352036.3 
LALA J142647.16+353612.6 
LALA J142630.34+354022.5 


5.746 
5.700 
5.674 


1.87 ±0.14 

3.88 ±0.15 
2.65 ±0.14 


11.7 

14 
8.8 


<360 

<460 
<220 


0.010 ±0.011 0.058 ±0.028 
0.013 ±0.012 0.115 ±0.029 
0.008 ± 0.012 0.018 ± 0.028 



^Units of 10 erg cm ^ s ^ . 

^Quoted uncertainties only account for random (photon counting) errors and exclude possible systematic 
errors. Red and blue side measurements are taken in 1200A wide windows beginning pa 30A from measured 
line center. 

Note. — Spectroscopic properties of the three confirmed z ?a 5.7 LALA sources. The redshift is based 
on the peak of the observed line profile, which may slightly overestimate the redshift since the blue wing of 
Lyman-a is absorbed. The measurement error on the line peak is estimated at Sz ^ ±0.002. The full width 
at half maximum (FWHM) is directly measured from the line profiles. The corresponding velocity width Av 
is based on subtracting in quadrature the effective instrumental resolution for a point source, and is therefore 
an upper limit since the galaxies may have angular size comparable to the 0.7" seeing prevalent in these 
data. The measured line fluxes agree very well with the photometrically inferred fluxes in table 1, except 
for LALA J142630. 34+354022. 5, whose line flux falls partially outside the narrow-band filter and is therefore 
underestimated from the narrow-band data. 
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Fig. 1. — Spectra of the three confirmed z 5.7 Lyman-a emitters. The top three panels show 
spectra smoothed with a Gaussian of 4A FWHM (comparable to the spectral resolution). The 
photon counting errors (accounting for the applied smoothing) are also plotted. Insets show the 
primary emission line, and are unsmoothed to allow line asymmetry to be studied. Finally, the 
fourth panel shows the night sky spectrum, reduced by a factor of 100 in order to plot it on the 
same scale as the data. (Night sky spectra for all three objects are in fact overplotted here; the 
differences among them are slight.) 
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Fig. 2. — Observational constraints on reionization. The solid circle shows the upper limit on the 
(spatially homogeneous) neutral fraction inferred from Lyman-a galaxy statistics at z ~ 5.7 (this 
paper and RMOl). The open squares with error bars show the evolution of the mass-averaged 
neutral fraction xhi that Fan et al (2002) derive from the Lyman-a forest and Gunn-Peterson 
trough in a sample of quasars up to redshift z = 6.28, using an inhomogeneous IGM model. 
Finally, the open circle shows the upper limit to the neutral fraction suggested by the discovery of 
a Lyman-ct-emitting galaxy at z = 6.56 (HCM 6a; Hu ct al 2002). Several loopholes might allow 
this source to be observed in a predominantly neutral IGM (e.g., Haiman 2002), but we regard 
these as unlikely (see section 4). 



